Critical behavior of an even offspringed branching and annihilating random walk 

cellular automaton with spatial disorder 
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A stochastic cellular automaton exhibiting parity conserving class transition has been investigated 
in the presence of quenched spatial disorder by large scale simulations. Numerical evidence has been 
found that weak disorder causes irrelevant perturbation for the universal behavior of the transition 
and the absorbing phase of this model. This opens up the possibility for experimental observation of 
the critical behavior of a nonequilibrium phase transition to absorbing state. For very strong disorder 
the model breaks up to blocks with exponential size distribution and continuously changing critical 
exponents are observed. For strong disorder the randomly distributed diffusion walls introduce 
another transition within the inactive phase of the model, in which residual particles survive the 
extinction. The critical dynamical behavior of this transition has been explored. 



I. INTRODUCTION 

The classification of the universality classes of nonequi- 
librium phase transitions to absorbing states is still an 
open problem of statistical physics [jj, |2j, |3j • Reaction- 
diffusion (RD) models exhibiting phase transitions to ab- 
sorbing states bear with a particular interest since many 
other types of systems like surface growth, spin systems 
or stochastic cellular automata can be mapped on them. 
Unfortunately there hasn't been experimental verifica- 
tion of such classes except the coagulating random walk: 
A + A — » A (CRW) in one dimension This is mainly 
due to the fact that the most well known, robust di- 
rected percolation (DP) class 5] is sensitive to disorder 
0, H CZ3| ! which occurs in real systems naturally. It 
would be very important to find some other nonequilib- 
rium class, which proves to be less sensitive to disorder, 
hence would provide a candidate for experimental verifi- 
cation. 

The study of disordered systems is a hot topic of cur- 
rent research of statistical physics ^lj ■ A principal con- 
dition for the relevancy of disorder is the Harris criterion 
|l2lll3| | set up for equilibrium systems and has been found 
to be valid in some nonequilibrium models. According to 
this criterion the pure critical point is stable against dis- 
order if the spatial correlation length critical exponent 
v±_ fulfills the inequality 



dv± > 2, 



(1) 



where d is the spatial dimensionality. However an excep- 
tion is reported very recently [l4| for DP with temporal 
disorder. Note that for CRW (which is exhibits the same 
scaling behavior as the A + A — > annihilating random 
walk (ARW) in Id) this criterion predicts relevant spa- 
tial disorder, (l>± = 1) still experiment p| did not report 
measurable effect of randomness unless very strong dis- 
order fractures the medium. 

Besides the robust DP an other well known universal- 
ity class is the so called "parity conserving" (PC) class 
of Id nonequilibrium transitions. This was discovered 
in a one-dimensional stochastic cellular automata (CA) 



exhibiting Zi symmetric absorbing states and domain 
walls following even offspringed branching and annihi- 
lating random walk: A — > 3A, 2A — > (BARW2) [l5|. 
Later it was observed by numerical studies of other mod- 
els d E El pi pi El El |H |H HI and field theo- 
retical studies [26l |27| confirmed the existence of a cor- 
responding fixed point distinct from that of DP. For a 
review see |2^| . This class is also called as directed Ising, 
DP2 or generalized voter model class. 

According to the Harris criterion disorder should be 
relevant for the critical behavior of this class = 
1.857(1) 3]). In contrast to this a recent renormaliza- 
tion group (RG) study Q did not find a strong disorder 
fixed point like in case of DP. The question naturally 
arises if BARW2 is really insensitive to disorder or the 
RG method 9] is not applicable for this case. The prin- 
cipal aim of the present study is to answer this question. 
Additionally in the absorbing phase of the BARW2 model 
the ARW dynamics dominates, which has also been ad- 
dressed in the studies j2j|[3(j ■ The renormalization study 
of ARW with spatial randomness in the reaction rates 
found margin al p erturbations to the fixed point of the 
pure system |30j. On the other hand an exact study 
of the infinite reaction rate ARW with space-dependent 
hopping rates found non-universal power-law decay of the 
density of A-s below a critical temperature [2^| . 

Note that in j^] the strong disorder is defined in such 
a way that it can not completely block reactions or dif- 
fusion of the reactants. Therefore the so called infinitely 
strong fixed point of does not correspond to the block- 
ing case. Such blocking or complete dilution was studied 
in a Id toy model of random quantum ferromagnetic Ising 
model |3l| where continuously variable power laws were 
found at the phase transition point. The effect of discon- 
nected domains in the reactions of CRW and ARW has 
been investigated in [32|]. This study reported stretched 
exponential decay in case of exponential domain size dis- 
tributions and continuously changing density decay for 
blocks distributed in a power-law manner. In the Id 
model we investigate such complete blocking may also 
occur, hence we investigate this topological effect. 
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II. THE NEKIMCA MODEL 

To study PC class transitions with disorder we have 
chosen a very simple stochastic cellular automaton (SCA) 
the NEKIMCA introduced in [2l|. It is easy to show 
that the dual variables of spins (f) the kinks (•) exhibit 
BARW2 dynamics via the synchronous spin-flip dynam- 
ics. In this SCA parity conserving kink branching is 
also generated due to the synchronous spin update of 
neighboring sites without introducing an explicit spin- 
exchange reaction as in case of the NEKIM model jlSj . 
The reactions are like the followings: 

• Random walk : f f • J, — ^ T • 1 I 

• Annihilation : | • I • j — ^ t T T 

w 2 

• Branching : | T • 4 I — ► T • 4 • T • I 

In the NEKIMCA there are two independent parameters 
parametrized as 

w % = T(l - S)/2 (2) 

w = r(i + g) (3) 

In the computer the state of a single spin is represented 
by a 1 or of a 32 or 64 bit word s(j) (depending on the 
CPU type). Hence 32 or 64 CA samples (exhibiting dif- 
ferent random initial conditions but the same quenched 
noise) updated at once. 

The following bit-parallel algorithm was used for the 
update of states s(j) at site j. A random number x(J) £ 
(0, 1) is selected with uniform distribution. If 

x(j) < qi(j) =Wi+x(j) (4) 

a spin-flip, corresponding to random walk of the dual 
variable 

s'(j) = [s(j + 1) A s ( j - 1)] A s(j) (5) 

is written to all bits of s(j)>. Following this another uni- 
formly distributed random number y(j) G (0, 1) is chosen 
and if 

y(J) < 9o(J) = w + xU) (6) 

a spin-flip, corresponding to annihilation of the dual vari- 
ables 

s'(j) = ([s(j - 1) A s(j)] k [s(j + 1) A s(j)]) A s(j) (7) 

is performed. Here x(j) denotes the quenched random 
noise variable with uniform distribution 

X(j)e(-£,e), (8) 

A and & are the logical XOR and AND of computer 
words. Note that for strong disorder qi{j) or q (j) may be 
less than at a site meaning a blockade for that reaction. 

A single Monte Carlo step (MCS) consists of updat- 
ing all s'(J) sites with periodic boundary conditions and 
writing back s(j) — s'(j) for j G (0, L— 1) (throughout 
the paper the time is measured by MCS). 
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FIG. 1: (Color online) Phase diagram of the disordered 
NEKIMCA for negative S. Bullets correspond to the dis- 
ordered PC class transition points, squares to the freezing 
transition (lines are for the guide of eye only). The dashed 
line shows w a , while dotted line denotes Wi. 



III. DYNAMICAL SIMULATIONS ALONG THE 
DISORDERED PC TRANSITION LINE 

The simulations were carried out on L = 4 x 10 4 — 10 5 
sized systems, with periodic boundary conditions. The 
initial states were randomly half filled lattices, and the 
density of kinks is followed up to 10 7 - 10 8 MCS. We 
started the exploration of the (5, e) phase diagram by 
determining the phase transition of the impurity free (e = 
0) case with r = 1. This was located at 5 C = —0.550(1) 
as a power-law decay of the kink density 

pozt- a (9) 

with exponent a = 0.28(1). This is in good agreement 
with the PC class value ElJ. Next we introduced the 
quenched disorder and determined S c for different e val- 
ues. Tableland Fig^Jshow the results. As one can see 



TABLE I: Numerical results for the disordered PC class tran- 
sition line. 



e 


Sc 


Q 





2 


V 


0.0 


0.550(1) 


0.280(6) 








0.1 


0.5513(5) 


0.280(6) 








0.2 


0.557(1) 


0.280(6) 








0.3 


0.5676(1) 


0.280(5) 


0.95(1) 


1.11(1) 


0.285(5) 


0.4 


0.5849(1) 


0.265(10) 








0.5 


0.6115(4) 


0.25(1) 


0.84(2) 


1.0(1) 


0.252(6) 


0.6767(2) 


0.7 


0.22(1) 


0.80(2) 







the disorder moves the transition point in such a way 
that the size of absorbing phase increases along the —S 
axis. 

Weak disorder (e < w ) seems to be irrelevant (see 
FigH}. The line corresponding to S — 0.5675 can well 
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FIG. 2: (Color online) Density decay of NEKIMCA 
for e = 0.3, T = 1 for different 8 
0.571, 0.57, 0.568,0.5675, 0.567, 0.0.565, 0.56, 0.55 (from 
top to bottom). The insert shows the local slopes for 
5 = 0.568, 0.5675, 0.567 (from top to bottom). 



be described by a power law from 10 3 < t < 10 8 MCS 
(5 decades !) with exponent a — 0.28. The local slopes 
(effective exponents) of the kink density decay are defined 
as 

-\n[p(t)/p(t/m)} 



eff(t) 



ln(m) 



(10) 



(where we used m = 2). At the critical point the a e //(t) 
curve exhibits a straight line shape for t — ► oo, while 
in sub(super)-critical cases a e ff(t) curves veer down(up) 
respectively. One can read off a e ff — * ol — 0.28(1) in 
a perfect agreement with the PC class value. Similar 
results are obtained for e = 0.1,0.2,0.4 and summarized 
in Table IJ 

To confirm these results we run cluster spreading sim- 
ulations from a single active seed (one kink, odd sector) 
and measured the number of particles (N(t)) and the di- 
ameter of clusters (R(t)). At the critical point these are 
expected to scale as 



N(t) oc f , R(t) oc t z ' 2 



(11) 



As Figure |21 shows both quantities exhibit power-law be- 
havior with PC class exponents: rj = 0.285(5), z = 

1.11(1) Qi. 

Finally we performed steady state simulations in the 
active phase at e = 0.3. We followed the density of kinks 
in many realizations until saturation is reached and aver- 
aged in a time window following that point. The steady 
state density in the active phase at a critical phase tran- 
sition is expected to scale as 

p{oo,8) oc \8-8 c f . (12) 

Using the local slopes method one can get a precise esti- 
mate for (3 as well as for the corrections to scaling 



Peff(§i) 



In p(oo, 8i) — In p(oo, Si 

Hh) - HSi-i) 



(13) 



FIG. 3: (Color online) Spreading in NEKIMCA for e = 0.3, 
r = 1 8 C = —0.5679. The insert shows the corresponding 
local slopes. 
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FIG. 4: (Color online) Effective order parameter exponent 
results on the disordered PC class line for e = 0.3, 0.5 (top to 
bottom). Lines correspond to linear fit. 



where we used the 8 C value determined before. The lo- 
cal slopes analysis resulted in /3 e // — > /3 = 0.95(1) in 
agreement with the PC class value again [25| (see FigQJ. 



A. Annihilation blocking on the disordered PC 
transition line 

In the inactive phase and at the critical point depletion 
of kinks dominates via the AA — > annihilation. How- 
ever for strong enough disorder active domains of arbi- 
trary sizes may also emerge due to the exponential dis- 
tribution of such events. The contribution of these such 
large regions can be estimated as in ljj with the differ- 
ence that for PC class even the absorbing phase decays 
algebraically i.e. one does not have exponential decay 
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that could slow down to stretched exponential at some 
"clean critical point" as in the DP case. The probability 
p a for finding a rare region of size l a is 

Pa oc exp(-ql a ) , (14) 

where q oc wf — w a + e is the probability that a site is 
active. The long-time kink decay density is dominated 
by these rare regions. For long times any finite region 
decays exponentially hence 

p{t) oc J dlj a p a exp(-i/<r(Z ) , (15) 

where r(Z a ) is the characteristic decay time of a size of re- 
gion of size l a . The average lifetime of such active regions 
grows as 

t(Z q ) oc exp(aZ Q ) (16) 

because a coordinated fluctuation of the entire region is 
required to take it to the absorbing state. The saddle 
point analysis of fTK)) as in results in power-law decay 

p{t) oc t- p "' a . (17) 

Hence one may expect continuously changing decay ex- 
ponents either as the effect of rare active regions or as 
the effect of high diffusion barriers [2^| . 

As we discussed in the previous section we did not 
see the occurrence of such continuous decay variables for 
weak disorder. However by going above e > ~ 0.4 the 
situation changes and we begin to see deviation from the 
pure PC decay behavior (see Table |]J The effect of dis- 
order becomes drastic. By increasing e above w the q 
reaction probabilities may become zero at certain sites 
randomly (with probability p w ) and the kink annihila- 
tion is blocked at those sites. From annihilation reaction 
point of view the systems falls apart to I blocks with an 
exponential probability distribution 

P(l) = PwPnw , (18) 

where the no-wall probability p nw is related to e accord- 
ing to our algorithm as p nw — 1 — p w = 1 — (e — w ), i.e 
where q < is set. This means that neighboring kinks 
can't even annihilate at the block boundary points which 
may cause relevant perturbation to the critical behavior. 

As Figure [5] shows for e = 0.5 the phase transition oc- 
curs at S c = —0.6115(3) and the density decay exponent 
is smaller than the PC class value: a = 0.25(1). For 
5 = -0.7, e c = 0.677(1), it is even smaller: a = 0.22(1) 
suggesting continuously changing exponents along this 
transition line by increasing e. Similarly the steady state 
exponent (3 and the spreading exponent z changes in this 
region (see Tableland Fig lT3|) . 



IV. DYNAMICAL SIMULATIONS IN THE 
INACTIVE PHASE 

In the inactive phase annihilating random walk dom- 
inates, in which for pure systems the density decays 
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FIG. 5: (Color online) Effective a in the strong disorder region 
of the disordered PC class transition line (e = 0.5), — S c = 
0.613, 0.612, 0.611, 0.609 (top to bottom). 
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FIG. 6: (Color online) Density decay in the inactive phase 
at coordinates: (8 = -0.6, e = 0.57), (8 = -0.5, e = 0.3), 
(8 = -0.525, e = 0.1) top to bottom). 

asymptotically as a power-law |33| 

pcct- 1/2 , (19) 

hence one can check the effect of disorder for this process. 
As one can see on Fig|S] the disorder does not change 
this behavior for small e, the ^(t)f 1 / 2 curves level off. 
For larger e however the deviations from the l|19(l law is 
observable for long times. Following a transient region, 
which is faster than 119|l and characteristic of ARW with 
finite reaction rates |2J] the decay slows down. For exam- 
ple for e = 0.57, 5 — —0.6 power-law fitting for t > 10 6 
MCS results in« = 0.43(1), while for e = 0.48, S = -0.5 
we obtained a = 0.46(1). These e dependent a values are 
in agreement with the results of |29j , where continuously 
changing power-law exponents are determined below a 
critical temperature in case of exponentially distributed 
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FIG. 7: (Color online) Density decay on the freezing transi- 
tion line for e = 0.6, 0.4, 0.2 (top to bottom). 

diffusion barrier heights. 



V. DYNAMICAL SIMULATIONS OF THE 
FREEZING TRANSITION 

As one can see on Figure ^ for negative 5 values the 
PC class transition line does not cross the Wi line mean- 
ing that the diffusion and branching reactions are not 
blocked. However due to the parallel computer algo- 
rithm we used in case of strong disorder (e > 1 — w ) 
q Q may become 1 at certain sites hence a random walk 
move is immediately overwritten by the J7J CA update at 
these sites. As the consequence in the absorbing phase 
(above the disordered PC transition line) the diffusion 
is blocked. Again the system falls apart to blocks with 
closed diffusion boundaries. As the depletion goes on 
kinks can not leave the blocks to interact with others in 
neighboring ones and "frozen" steady states emerge for 
e > —5. Note however that the system is not completely 
frozen, kinks can not only diffuse within the blocks, but 
branching may occur with (small probability) followed 
by a quick annihilation. In any case the residual density 
remains small. 

We run dynamical simulations on the freezing transi- 
tion line and as Figure0shows the density decay is slower 
than the square-root law of ARW l|19|) ) and changes con- 
tinuously along the freezing transition line (see Table ITT|) . 
However for e < 0.5 this deviation can be interpreted as 
a (logarithmic) correction to scaling. For e > 0.5 even 
the w a (annihilation) reactions are blocked and the de- 
viation from (|19|l is strong. We repeated this simulation 
on a system with different size (L = 10 4 ), but did not see 
any size dependence. 

The spreading exponent z changes similarly along the 
freezing transition line. It is close to z = 1 for e < 0.5 
and deviates it significantly for e > 0.5 (see Figure [SJ. 
The average number of kinks (N(t)) evolves to a small 



TABLE II: Numerical results for the freezing transition line. 



e 


-6c 


a 




z 


iV(oo) 


0.2 


0.2 


0.47(3) 


1.00(5) 


0.95(5) 


1.031(1) 


0.4 


0.4 


0.47(3) 




0.95(5) 


1.275(5) 


0.5 


0.5 


0.43(1) 




1.0(2) 


1.85(5) 


0.6 


0.6 


0.41(1) 


0.75(3) 


0.76(2) 


3.75(5) 


0.7 


0.7 


0.39(1) 


0.68(4) 


0.72(2) 


14(1) 
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FIG. 8: (Color online) Spreading in NEKIMCA along the 
freezing transition line e — 0.7, 0.6, 0.0, 0.2, 0.125 (top to bot- 
tom). The e = 0, 0.125 data correspond to T — 0.5 the others 
to r = 1. The e = 0.125 line corresponds to 5 — 0.5 the others 
to e = —5. The insert shows the corresponding local slopes. 

constant value as t — > oo (77 = 0) all along the freez- 
ing transition line. This value changes continuously by 
increasing the disorder strength. For e < 0.2 only a sin- 
gle kink survives (N(oo) ~ 1), while for other values see 
Table HH 

In the frozen (active) phase in the long time limit pairs 
within the same block annihilate and basically lonely 
kinks are wandering in confined regions of sizes I (lo- 
cal A —> 3A — > A reactions are also permitted with a 
small probability). Therefore the density of frozen kinks 
depend on the density of blocks, hence on e. The con- 
centration of blocks Cb can be expressed with the average 
length < I > of blocks as 

p oc c b oc< I > _1 . (20) 

For the exponential block size distribution 1|18|) the aver- 
age block size is 

<l>=(-\n(p nw ))- 1 (21) 
with the no-wall probability 

Pnw = 1 - Pw = 1 - (w + e - 1) = 2 - e - w . (22) 
The kink density is 

p oc -ln(2 - e-w ) , (23) 
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FIG. 9: (Color online) Effective order parameter exponent 
results on the freezing transition line for e = 0.2, 0.6 (top to 
bottom). Lines correspond to linear fit. 

which for small wall probability p w has the leading order 
singularity 

P oc p w ■ (24) 

This means that in the frozen phase = 1. The sim- 
ulations confirm this. For e < 0.5, where only diffusion 
traps are present the j3 ~ 1 indeed (see Table ITT1 and Fig- 
urcEJl. For e > 0.5 where annihilation is also blocked the 
exponent /3 decreases by increasing e (see Table CH. 

We have also investigated the phase space for positive 
S (with r = 0.5) where e > Wi causes a direct freezing 
transition via diffusion traps in the absorbing phase. The 
density decay simulations (see Figure ITU)) show small de- 
viations from the square-root law of ARW (|19fl ). which 
can be interpreted as corrections to scaling. 

VI. DISCUSSION 

We have investigated the phase diagram of the 
NEKfMCA as the function of quenched disorder with 
large scale simulations. We have chosen NEKIMCA, a 
very simple SCA, which is known to exhibit PC class 
transition of kinks to absorbing states. We determined 
static and dynamic exponents (a,[3,r],z) at several points 
along the phase transition line and in the inactive phase. 
The main consequence of this survey is that weak disor- 
der does not change the scaling behavior of the quantities 
studied. This supports the results of a recent real space 
RG study [9| but contradicts the Harris criterion. 

A possible resolution of this contradiction could be the 
study of the distribution of the critical point coordinates 
over the ensemble of samples of size L. According to 
recent progresses in the finite size scaling theory of dis- 
ordered systems, thermodynamic observable are not self- 
averaging at critical points when the disorder is relevant 
in the Harris criterion sense [H |H HHH HHH This 
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FIG. 10: (Color online) Density decay in NEKIMCA 
for positive S along the freezing transition line e = 
0.075,0.125,0.175,0.25 (top to bottom). The insert shows 
the corresponding local slopes. 



lack of self-averageness at criticality is directly related to 
the distribution of pseudo-critical temperatures T c (i,L) 
over the ensemble of samples (i) of size L. This has been 
shown numerically in case of a wetting transition of a 
polymer chain with quenched disorder |4l| . An inter- 
esting further direction of research would be the study 
of those distributions for nonequilibrium systems like for 
NEKIMCA. 

For weak disorder (e <~ 0.4) we did not see devia- 
tions from the asymptotic square root power-law (|19f) in 
the inactive phase. However for stronger disorder con- 
tinuously changing density decay exponent (a) were ob- 
served. This corroborates a former RG study [3(J for 
ARW with site dependent reaction rates and an exact 
calculation for the infinite reaction rate ARW with dis- 
ordered diffusion traps |29|. 

Very strong disorder in our model may introduce com- 
plete blocking of reactions or diffusion with exponential 
domain size distribution. If the system becomes seg- 
mented by diffusion walls a freezing transition to fluc- 
tuating states occur. In this case in odd parity blocks 
residual particles remain active. Large scale simulations 
suggest that this kind of disorder is marginal, the density 
decay slows down from the 1|19|) power-law by a logarith- 
mic correction or changes continuously if the annihilation 
is blockaded too. In the frozen phase the concentration 
of residual kinks increases logarithmically, with (3 e f / — >• 1 
asymptotically. This is in agreement with the results 
of [32|| , who considered the ARW model with complete 
blockades. However we did not see crossover to stretched 
exponential decay reported in for long times, but 
rather the decay slows down in our model. This is due 
to the fact that in our model in the frozen state not only 
blockades exist, but reaction and diffusion probabilities 
are randomized, which can introduce slower power-laws 
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In the region of the phase diagram where the disorder 
is so strong that the annihilation reaction becomes seg- 
mented we found continuously changing exponents along 
the disordered PC and the freezing transition lines both. 

Acknowledgements : 



We thank F. Igloi and G. M. Schiitz for the useful com- 
ments. Support from Hungarian research funds OTKA 
(Grant No. T046129) is acknowledged. The authors 
thank the access to the NIIFI Cluster-GPJD, LCG-GPJD 
and to the Supercomputer Center of Hungary. 



[1] For references see : J. Marro and R. Dickman, Nonequi- 
librium phase transitions in lattice models, Cambridge 
University Press, Cambridge, 1999. 

[2] For a review see : H. Hinrichsen, Adv. Phys. 49, 815 
(2000). 

[3] For a recent review see : G. Odor, Phase transition 
universality classes of classical, nonequilibrium systems, 
Rev. Mod. Phys. 76, 663 (2004). 

[4] R. Kroon, H. Fleurent and R. Sprik, Phys. Rev E 47, 
2462 (1993). 

[5] W. Kinzel, in Percolation Structures and Processes, ed. 
G. Deutscher, R. Zallen, and J. Adler, Ann. Isr. Phys. 
Soc. 5 (Hilger, Bristol, 1983); H. K. Janssen, Z. Phys. 
B 42, 151 (1981); P. Grassberger, Z. Phys. B 47, 365 
(1982); 

[6] P. Grassberger, Z. Phys. B 47, 365 (1982). 

[7] H. K. Janssen, Phys. Rev. E 55, 6253 (1997). 

[8] R. Cafiero, A. Gabrielli and M. A. Munoz, Phys. Rev. E 

57, 5060 (1998). 
[9] J. Hooyberghs, F. Igloi, Carlo Vanderzande, Phys. Rev. 

E 69, 066140 (2004). 
[10] T. Vojta and M. Y. Lee, Phys. Rev. E 72, 036126 (2005). 
[11] For a recent review see: F. Igloi, C. Monthus, Physics 

Reports 412, 277 (2005) 
[12] A. B. Harris, J. Phys. C 7, 1671 (1974). 
[13] A. J. Noest, Phys. Rev. Lett. 57, 90 (1986); A. J. Noest, 

Phys. Rev. E 38, 2715 (1988); 
[14] I. Jensen, J. Phys. A 38, 1441 (2005). 
[15] P. Grassberger, F. Krause and T. von der Twer, J. Phys. 

A:Math.Gen., 17, L105 (1984). 
[16] H. Takayasu and A. Yu. Tretyakov, Phys. Rev. Lett. 68, 

3060, (1992). 
[17] I. Jensen, Phys. Rev. E. 50, 3623 (1994). 
[18] N. Menyhard, J. Phys. A 27, 6139 (1994), 
[19] M. H. Kim and H. Park, Phys. Rev. Lett. 73, 2579, 

(1994). 

[20] D. Zhong and D. ben-Avraham, Phys. Lett. A 209, 333 



(1995) . 

[21] N. Menyhard and G. Odor, J. Phys. A. 28, 4505 (1995). 
[22] K. E. Bassler and D. A. Browne, Phys. Rev. Lett. 77, 
4094 (1996). 

[23] H. Hinrichsen, Phys. Rev. E 55, 219 (1997). 

[24] D. Zhong, D. ben-Avraham and M. A. Mufioz, Eur. Phys. 

J B 35, 505 (2003). 
[25] G. Odor and A. Szolnoki, Phys. Rev. E 71 066128 (2005). 
[26] J. L. Cardy and U. C. Tauber, Phys. Rev. Lett. 77, 4780 

(1996) ; J. L. Cardy and U. C. Tauber, J. Stat. Phys. 90, 
1 (1998) 

[27] O. Al Hammal, H. Chate, I. Dornic, M. A. Mufioz, Phys. 

Rev. Lett. 95, 100601 (2005). 
[28] N. Menyhard and G. Odor, Brazilian J. of Physics, 30, 

113 (2000). 

[29] G. M. Schiitz and K. Mussawisade, Phys. Rev. E 57, 
2563 (1998). 

[30] P. L. Doussal and C. Monthus, Phys. Rev. E 60, 1212 
(1999). 

[31] D. S. Fisher, Physica A 263, 222 (1999). 

[32] C. Mandache and D. ben-Avraham, J. Chem. Phys. 112, 

7735 (2000). 
[33] B. P. Lee, J. Phys. A 27, 2633 (1994). 
[34] L. Braunstein, H. O. Martin, M. D. Grynberg and H. E. 

Roman, J. Phys. A 25, L255 (1992). 
[35] S. Wiseman and E. Domany, Phys Rev E 52, 3469 (1995). 
[36] A. Aharony, A.B. Harris, Phys Rev Lett 77, 3700 (1996). 
[37] F. Pazmandi, R.T. Scalettar and G.T. Zimanyi, Phys. 

Rev. Lett. 79, 5130 (1997). 
[38] S. Wiseman and E. Domany, Phys. Rev. Lett. 81, 22 

(1998) ; Phys Rev E 58, 2938 (1998). 
[39] A. Aharony, A.B. Harris and S. Wiseman, Phys. Rev. 

Lett. 81, 252 (1998). 
[40] K. Bernardet, F. Pazmandi and G.G Batrouni, Phys. 

Rev. Lett. 80, 4477 (2000). 
[41] C. Monthus and T. Garel, J. Stat. Mech. (2005) P12011. 



